Crossover from critical orthogonal to critical unitary statistics at the Anderson 

transition 
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We report a novel scale-independent, Aharonov-Bohm flux controUed crossover from critical 
orthogonal to critical unitary statistics at the disorder induced metal insulator transition. Our 
numerical investigations show that at the critical point the level statistics are definitely distinct and 
determined by fundamental symmetries. The latter is similar to the behavior of the metallic phase 
known from random matrix theory. The Aharonov-Bohm flux dependent crossover is characteristic 
of the critical ensemble. 
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In analogy with classical phase transitions it is com- 
monly believed that the critical behavior at the disorder 
induced metal-insulator transition (MIT) is determined 
by fundamental symmetries. This is not only the case 
for the critical exponent of the localization length and 
the dc-conductivity, but also for the statistical proper- 
ties of energy eigenvalues and the wave functions. It was 
in particular the discovery of a new critical universal- 
ity class of spectral correlations ||,|^ that challenged our 
current understanding of the localization-delocalization 
transition. 

Random matrix theory (RMT) Q provides a gen- 
eral statistical description of the eigenvalue fluctuations 
and correlations in the metallic phase The corre- 
sponding universal energy level statistics only depend 
on fundamental symmetries of the system under consid- 
eration. For time reversal symmetry orthogonal (GOE) 
and symplectic (GSE) Gaussian ensembles are appropri- 
ate, whereas in case of broken time reversal symmetry 
the unitary (GUE) Gaussian ensemble describes the level 
statistics. In contrast to the metallic phase, the uncor- 
related energy levels of the localized states in the insu- 
lating regime are characterized by the Poisson statistics, 
independent of the symmetry. 

A critical statistics was not only found at the MIT 
in 3d in the presence of time reversal symmetry 0,^,^, 
but also for 2d symplectic systems [0J^ , where spin rota- 
tional invariance is broken but time reversal symmetry 
is conserved. In both cases the critical level statistics 
for small energy level separations closely resemble the 
distributions corresponding to the metallic phases. The 
probability density distribution Pc{s) of the normalized 
energy separation, s — \En~ En+i\/ A, is still Pds) 
for small s, where /? = 1 or 4 in the orthogonal and sym- 
plectic case, respectively. Here, i?„ and E^+i are two 
successive eigenenergies and A is the mean level spac- 
ing. This reflects the energy level repulsion due to the 
strong overlap of the corresponding eigenstates. From 
these results at the critical point and the universal be- 
havior in the metallic phase for all possible symmetries. 



one expects that Pc{s) ~ with /? = 2 for the critical 
unitary case, when time reversal symmetry is broken. 

Recently, however, it has been claimed |^ that the 
critical level statistics are independent of the presence or 
absence of time reversal symmetry. The authors based 
their assertion on numerical calculations of the nearest- 
neighbor energy spacing distribution Pc{s) and the A3- 
statistics (see below). This unexpected result seems to 
support an observation [ p^ that the critical exponent 
I' of the localization length is not changed by a mag- 
netic field. Because of the far reaching consequences for 
general theories concerning the MIT, a more careful and 
comprehensive study is indispensable. 

In this letter we present results of numerical investi- 
gations which demonstrate that there is a system size 
independent crossover between the critical orthogonal 
and the critical unitary statistics which is controlled 
by the magnitude of an applied Aharonov-Bohm (AB) 
flux. Besides the AB-flux model, we investigate two fur- 
ther mechanisms for breaking time reversal symmetry: 
a homogeneous magnetic field and a spatially randomly 
fluctuating magnetic flux with zero mean (random flux 
model). The same limiting critical unitary level spac- 
ing statistics P^{s) is found in all three cases. Our re- 
sults also show unambiguously that there exists a critical 
statistics for 3d unitary disordered systems that is clearly 
different from the critical orthogonal one. 

The dynamics of non-interacting electrons in a 3d dis- 
ordered system in the presence of processes that break 
the time reversal symmetry can be investigated by using 
the Anderson Hamiltonian 



H = J2e,\r){r\ + J2Vry\r){r'\ 



(1) 



The vectors r denote the sites of a simple cubic lattice 
with lattice constant a, and periodic boundary condi- 
tions are applied in all directions. The uncorrelated ran- 
dom energies are distributed with constant probability 
within the interval —W/2 < Cr < W/2, where W denotes 
the magnitude of the disorder. 
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The non-diagonal transfer matrix elements between 
nearest neighbors, 14, r' , contain the symmetry break- 
ing term: (i) for a constant magnetic field B, 14. r' — 
V exp{^i27ran), if r — r' = ±6^;, with r • ej, = n, and 
T4 r' — V else. The magnetic field is chosen to be com- 
mensurate with the lattice, (ii) for a random magnetic 
flux Vr.r' = y exp(i27r0r,r')i where the random phases 
0r,r' are drawn from an interval [—9o/2,9o/2] with uni- 
form probability p(^?r,r') — l/^'o- for the Aharonov- 
Bohm flux model 14, r' = V exp{±i27T<j>a/ L), r — r' — 
±ej, and j = x,y,z, where the AB-flux 0ab — 4>4>o is 
given in units of the flux quantum 0o = h/e, and L is 
the linear size of the system. V is the unit of energy. 
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FIG. 1. The critical unitary energy spacing distribution 
P"(s) in the presence of a magnetic field (a = 1/5) for sys- 
tems of size L/a — 5 (o), L/a = 10 (★), and L/a — 20 (o). 
For comparison, the critical orthogonal -P°(s) and the GUE 
result are also shown. The inset shows the large s behavior 
which can be fitted by an exponential decay exp( — ks) with 
K = 1.87 (dotted line). 

Energy eigenvalues were obtained for ensembles of sys- 
tems containing up to 32^ sites by direct diagonaliza- 
tion (using a Lanczos algorithm for large systems) or by 
means of the decimation method (l^Jl^ ■ The calculation 
was restricted to energy intervals about E/V = con- 
taining between 3% and 30% of the eigenvalues. Within 
this range we found no influence on our results. The 
number of realizations depends on the system size. The 
total number of eigenvalues amounts to 3 ■ 10^ {L/a — 5) 
and « 10^ for L/a = 10,20. A proper unfolding proce- 
dure has been applied to the energy spectrum in order to 
compensate for variations in the density of states. The 
normalized spacing distribution P{s) describes the fine- 
scale spectral structure and refiects for s ^ the level 
repulsion. We also calculated the Dyson-Mehta statistics 
A3 Q which measures the spectral rigidity and describes 
the spectral structure over a range of scales. 



In the presence of time reversal symmetry, the criti- 
cal disorder at {E/V = 0) corresponds to Wc/V ~ 16.4 
|^,|l3|,|5. It separates localized {W > Wc) from metal- 
lic states {W < Wc). We find that this critical disorder 
does not change when an AB-flux is applied. In contrast, 
Wc increases with increasing strength of a magnetic field 

We start to present our results with case (i), a 3d- 
Anderson model with applied strong constant magnetic 
field corresponding to a = eBa^/h — 1/5 fiux quanta 
per plaquette. The critical disorder is Wc/V — 18.1 
Ip^ . No dependence of the level spacing distribution on 
the system size could be detected within the error bars 
for L/a = 5, 10, 20 as expected for the critical ensemble. 
Similar results have also been obtained for magnetic field 
strengths corresponding to a = 1/4, 1/10 and 1/20. 

Fig. 1 displays the critical level spacing distribution 
P"(s) in the presence of a constant magnetic field in 
comparison with the critical data without magnetic field. 
The difference in the heights of the maxima and in the 
small-s behavior is clearly visible. This difference is 
much larger than the small fluctuations in the numer- 
ically obtained distributions. The small-s behavior is 
quadratic as expected for the critical unitary case, in 
contrast to the linear increase in the critical orthogonal 
ensemble. We tried to fit a suggested [0 critical level 
spacing distribution Pakl(s) = As'^ exp{—Bs'^) to our 
data with 7=1 + l/(i/d), where v and d are the critical 
exponent of the localization length and the Euclidean 
dimension, respectively. While it is possible to fit the 
bulk of our data with 7 w 1.2 which leads to v ^ 1.66, 
the deviation for larger s is obvious. The large-s behav- 
ior, which is shown on a logarithmic scale in the inset of 
Fig. 1, is best fitted by an exponential decay ~ exp(— ks), 
with K « 1.9 having approximately the same value as in 
the critical orthogonal case |]l^. A similar exponential 
decay with k « 4 has been reported for the 2d symplectic 
system 0. 

The Aa-statistics is defined as the least square devi- 
ation of the integrated density of states, N{E), from a 
linear behavior, averaged over an energy range k times 
the mean level spacing A centered at E 

A3(fc) = /mini j {N{E + e)-A-Bsfde\, (2) 

\ -fe/2 ' 

where ( ) denotes an ensemble average. The result shown 
in Fig. 2 for an energy interval around the critical point 
is again independent of the system sizes studied. The 
rigidity of the spectrum is enhanced in the critical uni- 
tary as compared to the critical orthogonal ensemble. 
The difference between the two critical curves is signif- 
icantly larger than the error bars, but smaller than the 
difference between GOE and GUE. The influence of the 
magnetic fleld to make the spectrum more rigid in the 
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metallic regime, is weakened at the critical point due 
to stronger disorder. Our results for P{s) and the A3- 
statistics clearly indicate the existence of a critical uni- 
tary statistics in the case of an applied strong magnetic 
field. 




FIG. 2. The Aa-statistics of the critical unitary ensemble 
(Wc/V = 18.1, Q = 0.2) for system sizes L/a = 5 (o), 
L/a = 10 (*), and L/a — 20 (o) in comparison with the 
critical orthogonal case, L/a — 20 (•). The known behavior 
of the insulating (Poisson) and metallic phases (GOE and 
GUE) are also shown. 



In order to investigate whether or not a weak magnetic 
field changes this behavior one must consider a different 
model, because due the requirement of commensurabil- 
ity the system size has to be taken to infinity as i? — > 0. 
Therefore, we used the random flux model (ii) which al- 
lows to study the effect of weak local magnetic fields with 
global magnetic field being zero. The strength of these 
local magnetic fields can be tuned by varying 9o, the 
width of the interval from which the random phases are 
chosen. We find the same critical unitary distribution 
P"(s) as in the strong magnetic field case, provided the 
condition L/lf, ;» 1 is fulfilled. The magnetic length 4 
which serves as a measure for the strength of the random 
local magnetic fields can be defined in analogy with the 
usual magnetic length = h/ {eB) — / (27ra). It is re- 
lated to the standard deviation, h, of the local magnetic 
fields. 



ll ^ h/ieb) = 



2^(l/i'Er(*r/<^0)' 



.1/2' 



(3) 



where the magnetic fiux <I>^ through a plaquette at site 
r, e.g. within the xy-plane, is calculated from the corre- 
sponding random phases $P = (/)o(^r+ex.r + ^^r+e^+ey.r + 

^r+e ,r+ex+e +^r,r+e )■ The total flux per cross-scction 



perpendicular to each of the three lattice axis is $t 



We now turn to the Aharonov-Bohm fiux model (iii) 
for which the absence of a distinct critical unitary statis- 
tics has been asserted |9). Here, we have to distinguish 
whether a magnetic flux is applied along one, two or 
three perpendicular directions. Depending on the magni- 
tude of the flux and the number of flux directions we ob- 
serve different scale independent level statistics for each 
set of parameters. To illustrate this unexpected behav- 
ior, the change of the maximum Pm{4') of the critical 
level spacing distributions is shown in Fig. 3, where the 
normafized maximum T^s,x{4>) = {PmW - Pm)/{Pm - 
P^) is plotted versus the Aharonov-Bohm flux (j). We 
find a smooth crossover starting from the maximum of 
the critical orthogonal P^ to the maximum of the criti- 
cal unitary ensemble P^ when the fiux is increased along 
two or three directions. Hence, in the limit of large 
flux there exists the same critical unitary level statis- 
tics f "(s) known from the random flux model and from 
the strong magnetic field case also for the 2- and 3- 
component AB-flux model. For one flux line in a sin- 
gle direction, however, this critical unitary distribution 
is not reached. Increasing the magnetic flux beyond 
(j) = 0.25 decreases Pm{4>) because of a ^o/2-periodicity 
in the flux: aX (f) — 0.5 all phases in Vry are equal to 
unity and the corresponding Hamiltonian is real. This 
behavior is known as "false T-breaking" ||l^. In all cases, 
within the error bars, the results are independent of the 
system sizes investigated. We find a similar crossover 
also for the As-statistics. 
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FIG. 3. The normalized maximum, rniax(0), of the critical 
level spacing distribution versus the Aharonov-Bohm flux (j) 
showing the system size independent crossover from critical 
orthogonal to the critical unitary ensemble. For AB-flux ap- 
plied in three perpendicular directions data are shown for 
L/a = 5 (O), L/a = 10 (+), L/a = 20 (□), and L/a = 32 
{•k). For AB-flux in two (o) and one (•) direction only the 
results for L/a — 5 are shown. 

The small-s behavior of the critical P(s, 0) for various 
values of the AB-flux applied in all three directions is 
shown in Fig. 4. While for small the overall deviation 
of the level spacing distribution from the critical orthog- 
onal curve is hardly visible, the small-s dependence is 
clearly quadratic. This holds for all fluxes considered, 
although the range of s values showing the s^-behavior 
decreases with decreasing flux. Only for (j) = the lin- 
ear relation, P{s) ^ s, characteristic of the orthogonal 
ensemble, is obtained. 
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FIG. 4. Double logarithmic plot of the small-s behavior 
of P{s) for various Aharonov-Bohm fluxes — 0.0 (•), 
(/) = 0.025 (*), cj> = 0.05 (o), (l> = 0.125 (□), = 0.25 (o). 

In the metallic phase, only the basic symmetry de- 
termines the characteristics of the spectral fluctuations 
even though the mechanisms to break the time reversal 
symmetry differ. In a magnetic field each closed electron- 
orbit not lying in a plane parallel to the magnetic field 
direction picks up a magnetic flux, which is in contrast 
to the smaller set of those closed trajectories with wind- 
ing number larger than zero, e.g. trajectories that circle 
round the single flux line in the AB-flux model. The 
latter situation was shown to be the reason for a flux 
controlled GOE — > GUE crossover in Aharonov-Bohm 
chaotic billiards |2^. Nevertheless, this crossover be- 
comes abrupt in the semiclassical limit L — > cx). A sim- 
ilar discontinuous transition is known from RMT 131 for 



infinite order matrices and also from numerical studies 
on metallic samples [ plf where the conductance diverges 
with L. 

Our results for the MIT show that the critical uni- 
tary statistics depend on how the time reversal symme- 
try is broken. This is due to the finite critical conduc- 
tance in combination with the system size independent 
Aharonov-Bohm flux 0ab = 4>h/e. Hence, the crossover 
becomes scale invariant at the metal-insulator transition 
point. This is different in the magnetic field case, be- 
cause here the magnetic flux (f>a = BL^ increases with 
system size. 

In conclusion, we have investigated the energy level 
statistics at the critical point of a 3d Anderson model 
when time reversal symmetry is broken by either a con- 
stant magnetic field, a spatially fluctuating magnetic 
flux, or an Aharanov-Bohm flux. A critical unitary level 
spacing distribution has been found that is distinct from 
the critical orthogonal one. We further showed that 
there exists a scale invariant AB-flux controlled crossover 
regime so that the critical level statistics depend on how 
the time reversal symmetry is broken. 
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